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Abstract 

In computing the third order terms of the series of powers of the center 
manifold at an equiUbrium point of a scalar delay differential equation, 
some problems occur at the term 10212:^2. More precisely, in order to deter- 
mine the values at 0, respectively — r of the function W2i{ . ), an algebraic 
system of equations must be solved. We show that the two equations are 
dependent, hence the system has an infinity of solutions. Then we show 
how we can overcome this lack of uniqueness and provide a formula for 

lf2l(0). 
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1 Introduction 

We consider the delay differential equation 

i{t)^ f{x{t),x{t-r)), (1) 

where the function / is a scalar function, /(O, 0) = 0, and r > 0. For the sake 
of simphcity we assume / S C°°{D) ((0,0) e C M^, I? open) but the results 
below are valid also for a lower class of smoothness (e.g. C'', k > 3). 

As usual in the study of delay differential equations [5], [3] we consider the 
function space B — : [— r, 0] 1— > R, is continuous on [— r, 0] } and its complex- 
ification, Be = B + iB. For a continuous function x : [— r, T] 1— ^ R, T > 0, we 
denote by Xt the element of B defined by Xt{s) — x{t + s), s e [— r, 0]. 

We also assume that the function / is such that, for every (j) £ B the above 
equation, with the initial condition xq = (j) has an unique solution x(t, cj)), t > 
and, thus, the equation generates a semi-dynamical system {T{t)}t>o on B, by 
T{t)^ = xt{^). ^ ^ 

In order to simplify the sequel developments, we write eq. ^ in the form 

x{t) = Ax{t) + Bx{t - r) + f{x{t), x{t - r)), (2) 



1 



/ being the nonlinear part of /, 

j,fc>0j + fc>2 

The hnearized equation attached to the previous, that is 

x{t) ^ Ax{t) + Bx{t- r), (3) 

generates a semigroup of operators on B. The eigenvalues of the infinitesimal 
generator of this semigroup are the solutions of 

\-A~ e-^'B = 0. (4) 

We assume that Q has a pair of pure imaginary complex conjugated solutions 
Ai_2 — ±wi, with w > 0, and all other eigenvalues have negative real part. 

Consider the subspace of Be spanned by (^1^2 and denote it by A4. The 
space Be is decomposed as a direct sum Be — Ai ^Af with the help of a 
projector. This projector is built with the help of a bilinear form, and this one 
is constructed by using the "adjoint" problem [1], [3]- 

In the above conditions, it is proved ^ that there is a local invariant mani- 
fold, called the center manifold that is a smooth manifold, tangent to the space 
Ai at the point x — and it is the graph of a function w(-) defined on a neigh- 
borhood of zero in A4 and taking values in Af. A point on the local invariant 
manifold has the form uipi +uip2 -\-w{uipi -\-uip2). 

In order to approximate the function w that defines the center manifold, we 
define w{z,'z) :— w{zipi +zip2). For simplicity, we will drop the hat, and we will 
denote also by w the new two-variables function defined above. 

We write ^ 

w{z,z)= —w^^jz'z^ (5) 

i+j>2 

where Wi^j S Be- From the invariance of the center manifold, a set of differential 
equations for the functions is obtained as well as a set of conditions for 
the determination of the integration constants. More precisely, these latter 
conditions have the form of linear relations between Wi,j{Q) and Wij{—r). 

For the second order terms of the series, the computation does not present 
any problems. This happens also for the third order terms, excepting u'2,1. For 
the function 1/72,1 the determinant of the linear algebraic system in Wij{0) and 
Wij{—r) is zero - Proposition 3.1. We prove that the system has (an infinity of) 
solutions - Proposition 3.2. 

The new problem is how to select the proper values for 102,1 (0) and u;2,i (—?'). 

In Section 4 we consider a perturbed problem, depending on a small pa- 
rameter e > 0, such that, when e — ^ 0, the unperturbed problem is obtained. 
We show that for this perturbed problem the system of algebraic equations 
for the corresponding Wc2,i(0), ^£2,1 (—7") has unique solution and that W2,i(0), 
W2,ii—r) can be obtained by taking the limit when e — >■ - Proposition 4.1 and 
Proposition 4.2. A consequence of the proof is the fact that the limit does not 
depend on the particular perturbed problem considered. A formula for W2,i(0) 
is given in Proposition 4.3. 
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2 The reduction of the problem to the center 
manifold 

For the presentation of our results we need the fohowing framework of [5] , 
[3], [I]- More precisely, we intend to describe the method for constructing the 
restriction of the problem to the center manifold. 

For this, we also need the space 

Bq — : [— r, 0] n- M, V' is continuous on [— r, 0) A 3 lim ij{s) G M > , 

and its complexification, Bqc- Notice that Bq consists of functions of the form 
ip ^ (fi + ado, where ip E B, a E M. and do : [—r, 0] i— )■ R, 



do{s) 



0, se[-r,0), 

1, s = 0, 



with norm given by ||?/;|| — \ip\Q + \a\. 

The linear part of the RHS of eq. ^ can be written with the help of linear 
operator L : B t-^ R, given by L{(p) = A(p{0) + Bif{~r). For further use, we 
write this operator, with the help of a Stieljes integral as 

J —r 

where 

r -B, s = -r; 
riis)^\ 0, .se(-r,0); (6) 

[ A, s = 0. 

In [U, the linear operator A : C^H-r, 0],M.) C Bq ^ Bq 

A{^) ^ ^ + do[L{p) ~ (7) 

is defined and it is proved that this is the infinitesimal generator of the semigroup 
of operators {S{t)}t>o given by S{t){(j)) = xt{(f>), where x{t,(j)) is the solution 
of equation ^ with the initial condition xq — (j). Then the nonlinear equation 
may be written as an equation in Bq, that is 

^=A{xt)+dof{xt), (8) 

where /((^p) = fifiO), v{~r)). 

In [3], the adjoint equation (associated to the linear equation (jS])) is 
defined as 

y{s) = -Ay{s) - By{s + r). 
The corresponding characteristic equation is 

X + A + Be^'' ^ 0, 

and it is obvious that, together with eq. it admits the solutions zLwi. The 
corresponding eigenfunctions are V'i(C) — e"'^*'', V'2(C) — C £ 
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Also in |2] , [3] in order to construct a projector on M , the following bilinear 
form is defined on C([0, r], C) x Be 



(V', ^) = ^(0)^(0) - / / V(C - 0MC)dCdv{9) = (9) 

J-r Jo 
J — r 

{rj being the function defined in ([5])). Then linear combinations of the functions 
ipj, j = 1,2, denoted by 5*^, i = 1,2 are constructed such that {'^i,ipj) = Sij. 
For this we determine the 2x2 matrix E, with elements e^j ~ {ipi, ipj): 



en = (V-i, V'l) = MO)fiW + B e-'^^'+^V"^" de = 1-{A- iuj)r, 

J — r 

ei2 = (V-i, = ^i(0)v52(0) +B f e-^'^^'+^^e-'^'de = 0, 
621 = (^2, V'l) = MO)MO) + B f e^'^^'+^h'^'de = 0, 

J ~r 

622 = (^-2, ¥^2) = V'2(0)¥'2(0) + S / 6^'^'-''+^'^ E-'^'^ dO = I - {A + lio)r. 

J —r 

Then 



^2 V ^^2 



and we obtain 



The projector defined in on Sqc £^nd with values in M. is given, for 
-0 = + doo- e Soc by 

T'W = ((^-l, + ^'l(0)(7) (^1 + ((^'2, 0) + ^'2(0)(7) (^2. (11) 

If V = e Sc, (cr = 0) we have = (l-i, 0)(^i + (^'2, 0)'^32- 

Now, for the solution z{-,(j)) of ([5]) with initial condition xq = (p, we can 
write 

Zt = '^lUl(t) + ip2U2{t) + V{t), 

with ui(t) = (^-i, Zt), U2(i) = (^'2, Zt), V = (/ - 7')zf 
In [1] the equation (|8l) is projected by T', to obtain, 

^ = Au+^'(0)/(<y9iUi+(^2U2+v), (12) 

dt 

where 

-^::)--(^o :).-(!: 
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iU w 



Here ui = U2, and the two scalar equations comprised in (jl2p are complex 
conjugated one to the other. Hence it sufhces to study one of the them. We 
denote ui = u and the projected equation is 

du ~~' 

— = {ti + iuj)u + ^i{0)f{ipiu + (p2U + v). (13) 
at 

If the initial condition is taken on the center manifold, then its image 
through the semigroup {T(i)}t>o, T{(j)) — Xtitp), is still on the manifold. Hence 

r(t)(0) = U{t)ipi + ^ip2 + W{u{t)ipi + W)V2), (14) 

with u(-), solution of the equation 
du — 

— = AlU + ^'l(0)/((y3lU + ip2U + w{uipi +ULP2)), (15) 

with the initial condition u(0) — uq, where V{4>) — uo^fi +?Io<^2- The real and 
the imaginary parts of this complex equation, represent the two-dimensional 
restricted to the center manifold problem. 

Now, if we use the series of powers of w from ([5]), we can write 

f{ipiu + ip2U + 'w{uipi+uip2)) ^ -T-rfij^ 
and, by putting 

5^,j=*i(0)A„ (16) 

eq. becomes 

du \ - 1 j 

— =Aiu+ 2^ —g,,,uv?. (17) 
Remark that = j , since the function / has real values. 

3 Computation of the coefficients j 

For the determination of the functions Wi.j of (O, the following relation 
(obtained from the invariance of the center manifold) is used [7] , [5] 

^ E -h^^^.M^'^^^ E -^5xfc"''^Vi(s)+ (18) 

j + k>2 ■' j+k>2 ■' 

j+k>2-'' ' j+k>2-'' 

This relation yields, by equating the terms of the same degree, differential equa- 
tions for each Wij . The integration constants are obtained from the following 
relation, also by equating the terms of same degree: 

Jt E -^"'^■■fe(0)"'"'+ E -^3.,fc"^^Vi(0)+ -^ff.,fc"'"V2(0) = 

i+k>2-'' ■ j+k>2-'' ■ j+k>2-'' 

(19) 

j+k>2-'' ■ j+k>2-'' ' j + fc>2-'' 
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3.1 The system for W2^i{—r), u;2,i(0) 

By matching the terms that contain u^u, we get the differential equation 
and the condition for W2^i. These are 

-^W2,iis) = u}iw2.i{s) + 52,ie"*" + 31 26""'"* + 2w2.o(s)5i,i + 

as ' ' 

(■5)52,0 + 2wi,i(s)gi 1 + U'o,2(s)5o,2> 

and 

wiu'2,i(0) + 2w2,o(0)5i,i + 1(0)52,0 + 2wi,i(0)5ii + two, 2(0)^0, 2 +52,1 +51, 2 = 

= Aw2,i(0) + Bw2^i{-r) + /2,i. 
From these we obtain the system of equations for W2^i{Q) and W2.i{~i')'- 



- e-'^'^W2,m + W2A~r) = -ff2,ire-"*'- + ^5i,2(e""^ " e"""^)- (20) 
-2gi,ie-'^''' f z«2,o(e)e""*'rf^-(ff2,o + 25i,i)e-"^'- f w^^m^-^^'dO- 



-{uji - A)w2,i{0) + Bw2,i{~r) = 52,1 + 5i,2 - hs + 251,1^2,0(0) + 

+ (52,0 + 25i_i)u;i,i(0) + 5o_2ti'o,2(0). (21) 
The matrix of this system is 



-(Loi~A) B 



and its determinant in IS. :— ui ~ A — Be that is equal to zero, since ui is a 
solution of the characteristic equation (|4]). Hence we obtained 

Proposition 3.1. The matrix of the system of algebraic linear equations 
for ?«2.i(0) o.iT'd W2^i{~r) has null determinant. 

In this situation the system has either an infinity of solutions or no solutions 
at all. The following Proposition solves this problem. 

Proposition 3.2. The equations of the system (|20p - (l21l) are dependent. 

Proof. We denote the right hand side of ([20)1 by Ri and the right hand side 
of (pij) by R2. We have to prove that 



BRi — R2. 

We prove this by showing that the following relations hold: 



B 



-92,ire 



52,1 + 51,2 - /2,i; 



(22) 



(i?l) 
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-2B5i,ie-"''- / W2,ome-^''d9 = 25i,iu;2,o(0); (i?2) 



(0)e-'^»d0 = (32,0 + 25, i)^«i,i(0); (i?3) 







(Rl). We first replace 52,1 with ^'i(0)/2,i, and ^ with *i(0)/2,i and (Rl) 
takes the form 



B 



-*i(0)/2,ire-"*'- + Avl/i(0)/2,i(e"^^ - e'^^ 
2uj 



= *l(0)/2,l+«'l(0)/2,l-/2,l. 



If /2,i = 0, then (-R1) is proved. If not, after dividing the relation by /2,i, 
we obtain 

- Be-'^"'*i(0)r + S^j^(e""' - e""''') = ^-1(0) + *i(0) - 1. (23) 
We prove that 

vS", (f)]i _ 
- Be-'^"^*i(0)r = *i(0) - 1, and B ' (e""' - e"'^*'^) = ^-1(0). (24) 

2a; 

We write (IMl)i as 

^1(0) + Be-"*'^«'i(0)r = 1 
and remark that it can be re-formulated as 

^-1(0) + B «'i(0)e-"'(''+'')e'^"'d6' = 1. 

J —r 

But this last relation is nothing else than 

= 1, 

that is true. Similarly, for (j24p 9. we write it as 

*i(0) + S— ^(e"*"" - e-'^n = 0, 

and rewrite this as 

_ ro _ 

^i{Q)+B / ^-1 (0)e'^^(^+'^)e"*^d6' = 0, 

J —r 

that is, actually, 

(*2,<^l) = 0, 



that is true. Since the two equalities of (j24p are proved, by adding them, we 
obtain (1^ . 
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Relations {R2)-{R4). If 51,1 0, then relation {R2) holds. We assume 
91,1 0- Similarly, for the two following relations, we assume 32,0 + '^Vi 1 7^ 
and 5o_2 7^ 0. 

We see that, after dividing with the assumed non-zero coefhcients, each of 
these relations has the form: 

-Be^'^"' Wj^k{d)e'"^''^d0 ^Wj,k{0) (25) 

J —r 

{3, k > 0, j + k = 2) that can be written as 

/■° 

Wj- fe(0) + B e-'^'^^+'''>Wj^k{0)de = 0. (26) 

J — r 

But relation (1^^ is in fact 

{i^2,Wj,k) = 0, 

which is true, since each of the functions Wj^k belong to the complementary of 
Ai. Hence all relations (i?2) - (i?4) are proved. 

Since we proved relations (i?l) - (i?4), by adding them, relation (|22p follows, 
and with it, the conclusion of our Proposition. □ 

The problem that occurs at this point is how to choose some " proper" values 
of ^2,1(0), U'2,i(— r) from the infinity of solutions of the system ([20| - (|2T]) . In 
the next section we give a solution. 



4 How to compute 

We consider a perturbation of our problem, of the form 

x{t) = A^x{t) + B^x{t -r) + f{x{t), x{t - r)), (27) 
where A^, B^ depend smoothly enough on e > 0, lim ~ A, lim B^ = B, and 

e\,0 e\,0 

are chosen such that, for small enough e, the linearized problem attached to ([77|) 
admits the eigenvalues Aei^2 — Me ± ^e*, with /ij > 0, while all other eigenvalues 
have negative real part. 

Obviously, from the construction it follows that lim — 0, lim = uj. An 

example of perturbed problem is given at the end of this subsection. 

From the construction, it follows that problem ^7} admits an unstable man- 
ifold, tangent to the space Me spanned by the two eigenfunctions <fei,2{s) = 
^{p.,_±iuje)s ^ corresponding to the two eigenvalues Aei.2. 

The unstable manifold is the graph of a function defined on and 
taking values in a subspace AC of Be, complementary to A^^. 

A procedure, similar to that presented in Section 2 for the non-perturbed 
problem, is developed in order to construct the restriction of the problem to the 
unstable manifold: 

-the adjoint equation and the eigenfunctions of its linearized, i.e. the func- 
tions V'£i(s) = e~^''' , ipe2{s) — e"'*"*, s G [—r,0], are considered; 

-the corresponding bilinear form, denoted also by ( • , • ) is constructed, and 
the functions ^^2 S C([0, r],C) are computed such that {'i'ej,ipek) = Sj^k] 
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these are _ 

l + {Xe-A,)r 

(1 — A^r + ^er)2 + Lu'^r^ 
the projector Ve : Be ^ Me, such that I — Ve : Be ^ J^e, is defined by 



We remark that, by the construction of the bihnear form, and by the defini- 
tion of the function whose graph is the unstable manifold, we have: 

(*ej,«^e)-0, J = 1,2. 

With the same type of reasonings (coming from [Ij), we find that the problem 
reduced to the unstable manifold is 

dv — 

— = AeW + ^'£i(0)/(t;(pei +Vipe2 + We(vipei +Vipe2))- 

As we did in Section 1 for w, we consider the function 

and (dropping the hat for simplicity of notations), we write: 

We{v,v)^ r^^w.^jv'v^ . (28) 

The coefficients Wdj are found by solving differential equations coming from 
relations similar to ((T5)) . (IT^ with Wj,k, gj,k, ■■■ replaced by w^j^k, dej.k, ■■■■ 
The equation for We2,i is 

^'^'Y^''^ = (2Ae + A,)u;e2.i(s) +5e2,ie^'" +5,1 2e^'" + 2we2,o(s)5ei.i + 
as ' 

+Wel,l{s)ge2,Q + 2w<:i,i(s)g,i_i + 
while the condition to determine the integration constant is 

(2Ae + Ae)u;e2,l(0)+2We2,o(0)gil+Wel,l(0)ge2,0 + 2w£l,l(0)g,ll+WeO,2(0)g£0,2+5e2, 1+5,1,2 

= ^w,2,i(0) + BeW,2A-r) + /e2,i. 

By integrating the differential equation above between —r and 0, we find the 
system of equations for We2,ii—'r), We2,i(0) : 

-e-(^^^+^')-w,2A0) + w,2A-r) = — V5.2,i(e-^^'"-e"(2^=+^')'-)- (29) 

At + Ae 

-^5.i,2(e-^^'- - e-(2^=+^=)'-) - 25.i.ie-(2^=+^=)'- j' u;,2.o(0)e~(2^^+^)^d0- 
-(5,2.0 + 25,i,i)e-(2^»+^')^ r w,,Ad)e-^''''-''^'d0- 

J — r 
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{A, - 2Ae - A,) U)e2,l(0) + B,W,2.i{-r) = .9^2,1 + 9el.2 - fe2.1+ (30) 
+ 2We2,o(0)gil + Wel,l(0)5e2,0 + 2Wel,l (0)g^i_i + We0,2 (0)5^o,2 ■ 

The matrix of the system is 

_g-(2A,+A,)r I 

^ A, ~ 2A, - A, B, 
with determinant 

Ae = -B,e-(2^=+^')'- - + 2Ae + K, (31) 

that is different of zero because otherwise the number 2Ac+A(: would be an eigen- 
value, with real part equal to Sfi,,, that contradicts the fact that all eigenvalues 
have real part < /ie. 

We denote the right-hand side of the two equations by R^2, respectively. 

Proposition 4.1 When e — ^ 0, the coefficients of system ((29|) - (p0| tend to 
the coefficients of system 

Proof. The assertion of the proposition is obvious for all the coefficients of 
the unknowns and for the terms in the right-hand sides, excepting the first term 
from the right hand sides of (gOl) and (gH). 

But also for these the conclusion comes easily 

e-^o Ae -I- Ae 

= - lim 5^2, le lim — -; r = -c/2,ie "*''r, 

^^o'' £^0 (Ae -I- A,)r- 

since Ag -I- Ae = 2/i(e) 0, when e — >• 0. □ 

Now, the natural idea is to solve system (P^I - ipU)) and to compute the limit 
of its solution when e — 0. As a matter of fact, it is enough to find z«e2,i(0), to 
compute its limit, and, then, 'W2,i{—r) will be found from one of the equations 
(EHl), (EH- By denoting the right hand sides of eqs. (EH), with R^i, R^2, 
respectively, we have 

/n\ B^R^i — i?e2 
We2,l(0) = • 

We know that both the numerator and the denominator of the above ex- 
pressions tend to when e — >■ 0, but the next Proposition shows that we can 
overcome this problem. 

Proposition 4.2 For any e > 0, B^R^i — R^2 and Ae can be written in the 
form 

BeRel — Re2 — Me^^ll^); 
Ae = /ie/l2(e). 
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where the functions hj, j = 1,2 have finite limit for e — >■ 0, and lim /12(e) / 0. 



e->0 



Proof. First we deal with the determinant A^. We have, by using the 
characteristic equation (Ae = + B^e~^^^), 

A, = -Bee-(2^^+^^)'' + B.e"^''^ + A, + A, = 
= B,e-^^^ (1 - e-2''^'')+2M. = B,e-^^^2n,r + + ..)+2f,, = 



2! 3! 



Hence the assertion concerning A^ holds, the expression of /i2 being obvious 
from the above relations. 

Moreover, we see that lini /i,2(e) = 2 (Be"'^*''r + l) = 2ruji - 2rA + 2, that 

can not be zero, since w 7^ 0. 

Now, in order to treat the term B^R^i — R^2, inspired by the proof of Propo- 
sition 3.2, we write 



BeRel — Re2 = Ei + E2 + E3 + £'4, 



where 



El 



Ae + A, 



■5e2,i(e' 



-Aer-_ -(2Ae + Ae)r 



-B, 



)-{9e2,l+9el,2-fe2,l), 



E2 = -2Be5ei,ie- 



-(2A.+Ae)r 



W. 



'e2,o(^)e-('^'+^')^de - 23,i,iw;e2,o(0), 



E3 = -i3,(ffe2,o+25ei,i)e"^'^^+^^''' / i«ei,iWe-('^=+^^)'d0-(5.2,o+25,i,iKi,i(O), 



£^4 = -5^5.0,26-^'^'+^'^'" f «;.o,2We-^'^^+^'^'rf^-5.o,2«'.o,2(0). 

El. We write £1 as Ei = fe2,i{Eii + -E12), where: 
-B, 



E 
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■4',i(0)(e-^''' - e-(2^«+^«)'-) - *,i(0) + 1, 



E12 = — ^*,i(0)(e-^''- - e-(2^'+^')'-) - *ei(0). 



2A 



For £11, we have 



= 1- (*,l,77e), 

where ?7e(s) = e(2>*«+^«)*, s g [-r, 0]. 
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But we know that 1 = {^ei,ipei), and, substituting this above, we find 

^11 = {^el,Vel) - i^eUVe) = el, Vel - Ve) ■ 

But 



:i(s) - r]e{s) = e^^' - e(2^'+^'> = e^'^ (l - e^"'^) = (32) 



( 1 



2/i(;S (2/igS 



|2 



2! 3! 



The series in the paranthesis is convergent and its sum is a bounded function 
on [— r, 0]. We set 

P.(a) = -2.e-(l + ^+(^ + ... 

and we obtain 

Now we pass to £^12 



= -vl/,i(0) - fi,M/,i(0)e-^«'^ / e^^»^ds = 

J— r 

_ _ /-o - 

But (^'e2,¥'ei) = 0) and, then, we can write 

£^12 = {^e2,Vel) - (*e2,?7e) = e2, Vel - Ve) , 

and from this point, by repeating step by step the reasonings made for En we 
obtain 

£^12 = /Xe(*62, Pe)- 

Finally, 

El = Pefe2,l ((*el, Pe) + (*e2, Pe)) ■ (33) 

E2 — E4. For a unitary writing, we define 

"2,0 := -2.9el,l, :== -ge2,0 - '^Vel,!, "0,2 := -5e0,2- 

Then, each of the expressions Ei, i = 2,3,4, can be written as 

aj,k (^Wej,k{0) + y e-(2^'+^''(''+'')u;„-,fe(s)ds) = aj,k{Ve, Wej,k), 

where j, fc > 0, j + fc = 2, and rje{C) = e-(2->'=+^-)^ , C e [0, r]. 

From the definitions of the function We, that defines the invariant manifold, 
and of the projector on A4e, we have 

{lpel,Wej,k) = 0. 
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Then we may write 

and _ 

5.(0 - *,.(0 = e-<"-*'-K - e-».< = e-*.< - l] = 

We consider the function p : [0, r] i-> K, 

and, finally we can write B^R^i — R^2 — fJ'ehi{e), where 

hl{e) = /e2,l((«'.l, P£) + (*e2, Pe))-25rel,l(Pe,We2,o)- 
-(5£2,0 + 2g,i i)(p,,-U;,i,i) -5,o_2(Pe:W'eO,2)- 

The proof is complete. □ 

Now we can easily compute u'2.i(0)- 
Proposition 4.3. The value of W2,i(0) is 

/2.i(*i + *2, p) - 2gi,i{p,W2.o) - (.92,0 + 2gi i){p,wi,i) -902(9^^^,2) 
W2 i(0) = : ■ ■ . 

(34) 

where p{s) = -2se'^'\ s G [-r,0] and p{C) = -2(6""*^, ( e [0,r]. 
Proof. The preceding Proposition implies 

lim W(:2,i(0) = lim 7^4~T- 

In order to see that the coefficients g^j.k and the functions w^j^k in ^i(e) do not 
present any problem when passing to limit, we list their values in the Appendix. 

We see that limc^^o 5ej,fc = 5j,fc, and that w^j^k — >■ Wj.k when e — >■ 0, uniformly 
on [— r, 0]. Then, by observing that 

lim«,(s) = -2se'^^^ 

e->0 

the convergence being uniform with respect to s e [— r, 0] and 

limp,(C) = -2Ce-'^<, 

c— >-0 

(uniform convergence on [0,r]), the result of our Proposition is obtained. □ 

From any of the two equations of system (pHl) - (PT|) we then find 'W2,i{-~r). 
From the proof of Proposition 4.2, we see that lim 'Wc2 i(0) does not depend 

e-i-O ' 

on the specific perturbation of the problem chosen. 

Example of a concrete perturbation. We take = B(l + e) and chose 
in order to have two eigenvalues of the form /i(e) ± wi, with /i(e) > 0. 



13 



The characteristic equation associated to the hnear part of the equation is 



Taking the imaginary parts of the equation, we find 

B sin(a;r) 



In 



(35) 



(36) 



Since Lui is an eigenvalue for the hnearized of our problem, the characteristic 
equation for e = implies 

_Bsin(cLir) 

B sin(ajr) , , 
and, thus, for e > 0, ^ — -(1 + e) > 1 and ^(e) > 0. 



By taking the real part of we obtain 

= ^(e) - B{1 + e)e-^(") cos(a;r). 



(37) 



Problem (j27p with A^, B^^ given above is an example of perturbed problem. 
For small enough e, its linearized part has the eigenvalues Xei.2 = M(e) i ito, 
with /i(e) > 0, while all other eigenvalues have negative real part. 

Remarks. 1. In [7], for a particular problem, when trying to compute W2,i 
(denoted otherwise there), a perturbed problem (depending on a small param- 
eter e) is considered, the corresponding We2,i is computed and the limit when 
e — is taken. But, besides the fact that it is a particular problem (hence the 
result is not general), there is no proof there that the limit does not depend on 
the specific perturbation chosen. 



2. In a previous work [S], in which we intended to present an example of 
Bautin type bifurcation in a delay differential equation, we encountered the lack 
of uniqueness in determining W2,i(0), W2,i{—r). The equation considered there, 
that did not come from some model (it was constructed by us) was 

X ^ ax{t - r) + x^{t) + cx{t)x{t - r), (38) 

with r = 7r/2. Obviously, a; = is an equilibrium point, and we looked for the 
values of the parameters a, c where the sufficient conditions for Bautin type 
bifurcation are fulfilled [1]. We showed that the linear problem associated has, 
for ao — "1 two pure imaginary eigenvalues, Ai,2 = ±i, and that for two values 
of c, 

_ 18 - Ttt ± V36 + 2127r + 7r2 
" 2(37r-2) ' 

the first Lyapunov coefficient is zero. Up to this point the values of W2,i were 
not necessary. 

Then, we intended to compute the second Lyapunov coefficient since, as 
we remarked in [4], when this is positive, the two limit cycles (one inside the 
other) that occur in the Bautin type bifurcation, exist when two eigenvalues with 
positive real part exist, i.e. the cycles actually exist on the unstable manifold. 

When computing the second Lyapunov coefficient l2{c), we needed the val- 
ues of ^2,1(0), if2,i(— r), but since the two equations that yield these values are 
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dependent, at that moment we have chosen arbitrarily 102.1(0) = and com- 
puted W2,i{—r) from one of the two equations. We found that hici) > 0, while 
^2(02) < 0. 

However, the correct way of solving the problem is the one presented in this 
paper, and this Remark is intended to play the role of an Erratum to [5]. 

With the method developed here, by using formula we find, in the case of 
ci(w 1.52799), u;e2,i(0) = 0.285 - 0.28i, w;,24(-r) = 1.442 - 1.612i, h = 4.528. 
Hence, in the parameters plane, in a neighborhood of the point a = — 1, c = ci, 
there is a zone where the unstable manifold exists and for parameters a, c in 
a subset of this zone, two periodic orbits (one inside the other) exist on the 
unstable manifold. 

We re-analyzed the case of C2(~ —2.06554) and found: 
We2,i(0) = -0-69 - 0.278i, We2,i(-0 = -4.732 - 1.537i, h = 3.726. 

This shows that equation psp presents Bautin type bifurcation for both 
values ci, C2. 

5 Appendix 

First we write down the expressions of f^ij, i+ j = 2, and that of /e2,i- We 
have: 

/ei,i - C2.0 + Ci,i(e-^''- -I- e-^'"^) + Co.2e~^^'\ 

fe0,2 — /e2,0 = ^2,0 + 2Ci_ie ^'^ + ^0,26 
/e2,l = C2.O (2«;el,l(0) + W,2..0m + 

+Ci,i (w,2fime~^'^ + 2u;,i,i(G)e-^^'^ + zi'e2,o(-r) + 2w;,i,i(-r)) + 
+C0.2 (2w,i,,{~r)e-^^^ + w,2.A-r)e-^''') + ^3,0 + C2.1 (2e-^'- + e"^^'') + 

By multiplying the above quantities by 4'ei(0) we obtain the corresponding gei,j- 
Now we compute w^j^k, j + k = 2. The differential equation for We2.o is 

^^^^^^ = 2A,«;,2,o(s) +ff.2,oe^=^ +5.o,2e^'^ 
and by integrating it, we obtain 

^^.2,0(5) = w,2.o{Oy^'' ^—-9e2,o {e^'' - e^^'') - 



\.s 2\cS 

By taking s — ~r we get 



"'^^'^w^.2,o(0) + ^^e2.o(-r) = ^5.2,0 (e"^^'' - e- 
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while the supplementary condition is 

(Ae - 2Ae)w<:2,o(0) + BeWe2,o{-r) = 5^2,0 + 5e0,2 " fe2,0- 

The differential equation for is 
^^^^^ = 2Me^i;.i,i(s) 

from where, 

We set again s = —r, 

and by using also the condition for Wei,i : 

we have the two equations for Wci^i{—r), ■u;<;i,i(0). 
For Weo,2, it is enough to remark that We0,2 = We2,o- 
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